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Abstract 

In this paper we introduce a stochastic integral with respect to the solution 
X of the fractional heat equation on [0,1], interpreted as a divergence operator. 
This allows to use the techniques of the Malliavin calculus in order to establish an 
Ito-type formula for the process X. 
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1 Introduction 

In the last past years, a great amount of effort has been devoted to a proper definition of 
stochastic PDEs driven by a general noise. For instance, the case of stochastic heat and 
wave equations in M n driven by a Brownian motion in time, with some mild conditions 
on its spatial covariance, has been considered e.g. in jUEJE], leading to some optimal 
results. More recently, the case of SPDEs driven by a fractional Brownian motion has 
been analyzed in [3J m the linear case, or in EI] for the non-linear situation. 

In this context, it seems natural to investigate the basic properties (Holderianity, be- 
havior of the density, invariant measures, numerical approximations, etc) of these objects. 

♦Partially supported by the CONACyT grant 45684-F 
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And indeed, in case of an equation driven by a Brownian motion, a lot of effort has been 
made in this direction (let us cite ^3 E3 IE] among others). On the other hand, results 
concerning SPDEs driven by a fractional Brownian motion are rather scarce (see however 
[TT] for a result on SPDEs with irregular coefficients, and ^3] for a study of the Holder 
regularity of solutions). 

This article proposes then to go further into the study of processes defined by fractional 
PDEs, and we will establish a Ito-type formula for a random function X on [0, T] x [0, 1] 
defined as the solution to the heat equation with an additive fractional noise. More 
specifically, we will consider X as the solution to the following equation: 

d t X(t,x) = AX(t,x) + B(dt,dx), (t, x) e [0, T] x [0, 1], (1) 

with Dirichlet boundary conditions and null initial condition. In equation (pQ), the driving 
noise B will be considered as a fractional Brownian motion in time, with Hurst parameter 
H > 1/2, and as a white noise in space (notice that some more general correlations in 
space could have been considered, as well as the case 1/3 < H < 1/2, but we have 
restrained ourselves to this simple situation for sake of conciseness). 

Then, for X solution to JTJ, t e [0,7], x e [0, 1] and a Cf-function / : R -> R, we will 
prove that f(X(t,x)) can be decomposed into: 

f(X(t,x)) = f(0)+ f [ {MU'(X))(s,y)W(ds,dy) + l f f'{X{s,x))K x {ds), (2) 
Jo Jo z Jo 

where in the last formula, M£ x is an operator based on the heat kernel Gt on [0, 1] and the 
covariance function of B, W is a space-time white noise, and K x is the function defined 
on [0,7] by: 

PS PS 

K x {s) = H(2H - 1) / / G 2s - Vl - V2 (x,x)\v 1 - v 2 \ 2H ~ 2 dv l dv 2 . 
Jo Jo 

Notice also that, in (J2J), the stochastic integral has to be interpreted in the Skorohod sense 
(see Theorem 13 . 1 31 for a precise statement). 

As mentioned above, once the existence and uniqueness of the solution to (JTJ) is es- 
tablished, it certainly seems to be a natural question to ask whether an Ito-type formula 
is available for the process we have produced. Furthermore, this kind of result can also 
yield a better understanding of some properties of the process itself, such as the distri- 
bution of hitting times, as shown in j^J. It is also worth mentioning at this point that 
formula (0) will be obtained thanks to some Gaussian tools inspired by the case of the 
fractional Brownian motion itself. This is due to the fact that X can be represented by 
the convolution ^ 

X{t,x)= I [ M t>s (x,y)W(ds,dy) (3) 
Jo Jo 

of a certain kernel M on [0,7] x [0, 1], defined at (|19jh with respect to W. This kind of 
property has already been exploited in |Hj for the case of the heat equation driven by a 
space-time white noise, but let us stress here two differences with respect to this latter 
reference: 
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On the one hand, an important step of our computations will be to obtain the 
representation (jHJ) itself (see Corollary 13. 3j) and to give some reasonable bounds on 
the kernel M and its derivatives. 

On the other hand, the little gain in regularity we have in the current situation with 
respect to jHj will allow us to obtain a formula for t i— > f(X(t, x)), while in the latter 
reference, we had to restrict ourselves to a change of variable formula for 

1 1-> I f(X(t,x))i/i(x)dx, 



o 



for a continuous function ijj. 



Let us say now a few words about the method we have used in order to get our 
result: as mentioned above, the first step in our approach consists in establishing the 
representation (J3"|) . This representation, together with the properties of the kernel M, 
suggest that the differential of X should be of the form 



X(dt,x) 



o jo 



d t M t Jx,y)W{ds,dy) 



dt. (4) 



This formula is of course ill-defined, since (s,y) t—> d t M tjS (x,y) is not a L 2 -function on 
[0, t] x [0, 1], but it holds true for a regularization M e of M. We will then obtain easily 
an Ito type formula for the process X £ corresponding to M £ , where the differential P|) 
appears. Therefore, the main step in our calculations will be to study the limit of the 
regularized Ito formula when e — > 0. Notice that this approach is quite different (and 
from our point of view more intuitive) from the one adopted in PJ Ej , where the quantity 
E[f(X(t,x))I n (ip)] was evaluated for an arbitrary multiple integral I n {f) with respect to 
W. 

Our paper is divided as follows: at Section 2, we will describe precisely the noise and 
the equation under consideration, and we will give some basic properties of the process X. 
Section 3 is devoted to the derivation of our Ito-type formula: at Section 3.1 we obtain the 
representation (j3J) for X, the regularized formula is given at Section 3.2, and eventually 
the limiting procedure is carried out at Sections 3.3 and 3.4. In the sequel of the paper, 
c will designate a positive constant whose exact value can change from line to line. 



2 Preliminary definitions 

In this section we introduce the framework that will be used in this paper: we will define 
precisely the noise which will be considered, then give a brief review of some Malliavin 
calculus tools, and eventually introduce the fractional heat equation. 

2.1 Noise under consideration 

Throughout the article, we will consider a complete probability space (Q, P) on which 
we define a noise that will be a fractional Brownian motion with Hurst parameter H > 1/2 
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in time, and a Brownian motion in space. More specifically, we define a zero mean 
Gaussian field B = {B(s, x) : s G [0, T], x G [0, 1]} of the form 

B(t,x)= [ [ K H (t,s)W(ds,dy). (5) 
Jo Jo 

Here W is a two-parameter Wiener process and Kh is the kernel of the fractional Brownian 
motion (fBm) with Hurst parameter H G (|, 1). Namely, for < s < £ < T, we have 

s) = C H s^' H J (u - s) H ~iu H ~^du, 

where Ch is a constant whose exact value is not important for our aim. Observe that the 
standard theory of martingale measures introduced in easily yields the existence of 
the integral (jHJ). 

Note that it is natural to interpret the left-hand side of © as the stochastic integral 

B{l MxM ) ■= [ [ B(ds,dy). (6) 
Jo Jo 

The domain of this Wiener integral is then extended as follows: let TC be the Hilbert space 
defined as the completion of the step functions with respect to the inner product 

(1[o,«]j l[o,t])w = (KH(t,-),K H (s,-)) L 2( [0!T]) 

= H(2H - 1) f I \u- r\ 2H ~ 2 dudr. (7) 
Jo Jo 

Thus, by Alos and Nualart j2] , the kernel Kh allows to construct an isometry K* H T from 
WxL 2 ([0,l]) (denoted by H T for short) into L 2 ([0, T] x [0, 1]) such that, for < s < t < T, 

(K HtT l[ ,t]xlo,x]) (s,y) = K H (t,s)l[ 0tX ](y) 

= l[0p\(y) / 1 [o,t](r)d r K H (r,s)dr. 



Therefore the Wiener integral © can be extended into an isometry (p i— > B(ip) from T~Ct 
into a subspace of L 2 (fl) so that, for any tp G Ht, 

B(<p) = [ T f\K HtT cp){s,y)W{ds,dy). (8) 



o J o 



Then, for two elements (p and i/j of Ht, the covariance between B(<p) and -B(V') is given 
by 



E[B(<p)B(rl>)] = H(2H-l) [ [ 

Jo Jo 



T r T r l 

2H-2 



<p(s, y)\s — r\ ip(r,y) dsdrdy. (9) 



Notice that an element of Tix could possibly not be a function. Hence, as the in 
fBm case, we will deal with the Banach space \Ht\ of all the measurable functions ip : 
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[0, T] x [0, 1] -> R such that 

|M|| Wt , = H{2H-l)f [ [ \v(r,y)\\u-r\ 2H - 2 \v(u,y)\dydudr 

Jo Jo Jo 

= So So \Ss l^( r '^l^ r ^( r ' s ) rfr ) dsdy <oo. 
It is then easy to see that L 2 ([0, T\ x [0, 1]) C \H T \ C ft T - 

2.2 Malliavin calculus tools 

The goal of this section is to recall the basic definitions of the Malliavin calculus which 
will allow us to define the divergence operator with respect to W . For a more detailed 
presentation, we recommend Nualart [10J. 

Let S be the family of all smooth functionals of the form 

F = f(W(s 1 ,y 1 ),...,W(s n ,y n )), with (s hyi ) G [0,T] x [0,1], 

where / G C£°(R n ) (i.e., / and all its partial derivatives are bounded). The derivative of 
this kind of smooth functional is the L 2 ([0, T] x [0, l])-valued random variable 

i=i 1 

It is then well-known that D is a closeable operator from L 2 (Q) into L 2 (Q x [0,T] x 
[0,1]). Henceforth, to simplify the notation, we also denote its closed extension by D. 
Consequently D has an adjoint 5, which is also a closed operator, characterized via the 
duality relation 

E (F5(u)) = E ({DF,u) L 2 {[0 ,t M o,i])) , 

with F G S and u G Dom(<5) C L 2 (VL x [0, T] x [0, 1]). The operator 8 has been considered 
as a stochastic integral because it is an extension of the ltd integral with respect to W 
that allows us to integrate anticipating processes (see, for instance, [TU]). According to 
this fact, we will sometimes use the notational convention 

8(u) = / u s , y W(ds, dy). 
Jo Jo 

Notice that the operator 5 (or Skorohod integral) has the following property: Suppose 
that F is a random variable in Dom(D) and that u is Skorohod integrable (i.e., u G 
Dom(5)), such that E(F 2 J Q T J^(u(s,y)) 2 dyds) < oo. Then 



T rl rT pi pT pi 

Fu(s, y)W(ds, dy) — F \ j u(s,y)W(ds,dy) - / (D s>y F)u(s,y)dyds, 

Jo Jo Jo Jo 

(10) 



o </o 



in the sense that (Fu) G Dom(<5) if and only if the right-hand side is in L (f2) 
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2.3 Heat equation 



This paper is concerned with the solution X to the following stochastic heat equation on 
[0, 1], with Dirichlet boundary conditions and null initial condition: 



d t X(t, x) = AX(t, x) + B(dt, dx), (t, x) e [0, T] x [0, 1] 
X{0,x) = 0, X(t,0) =X(t,l) = 0. 



It is well-known (see [Ug\) that equation (|11|) has a unique solution, which is given explicitly 
by 



X(t,x) = / G t ^ s {x,y)B(ds,dy), 
Jo Jo 



where 



exp 



(y — x — 2n)' 
it 



exp 



(y + x — 2n)^ 
4t 



(12) 



(13) 



stands for the Dirichlet heat kernel on [0, 1] with Dirichlet boundary conditions. Let us 
recall here some elementary but useful identities for the heat kernel G: 

Lemma 2.1. The following relations hold true for the heat kernel G given by A 131) : 



G t (x, y)dy = 1, G t (x, y) < exp 



c 2 (x - yf 



and 



\d t G t {x,y)\ < -^exp 



Cj{x - yf 
t 



for some positive constants Ci, C2, C3 and C4. Furthermore, G can be decomposed into 

G t (x,y)=G 1>t (x,y) + R t {x,y), (14) 

where 



Gi, t (x,y) 



/Ant 



exp 



(y - xf 

At 



exp 



(y + xf 
At 



— exp 



{y + x-2f 
At 



and R t (x,y) is a smooth bounded function on [0, T] x [0, l] 2 . 

Let us recall now some basic properties of the process X defined by (JTTJ) and (jEJ), 
starting with its integrability. 

Lemma 2.2. The process defined on [0, T] x [0, 1] by hl2i) satisfies 

sup 

te[o,r|,xe[o,i] 
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Proof. We have, according to Q and Lemma \2. 11 that 

E[\X(t,x)\ 2 ] = c H i2=2H / G t-s(x,y)G t - u (x,y)dy 

J\0,t] 2 \ s ~~ u \ Jo 



'[<>■']■-' 

< c / 7- r^-r / G t -u(x,y)dy 



s — u 



c 



dsdu 

0,t] a (*-s)V2| s _ M |2 



-2H : 



and the last integral is finite by elementary arguments. □ 

One can go further in the study of X, and show the following regularity result (see 
also [HJ): 

Proposition 2.3. Let X be the solution to Ml\) . Then, forti,t 2 G [0,T] and a; G [0,1], 
we have 

E[\X{t 2 ,x)-X(t 1 ,x)\ 2 ] <c\t 2 -t x \ 2 \ 

for any 7 < H — 1/4. In particular, for any T > and a; G [0, 1], tae function t G [0, T] 1— > 
X(t, x) is 7 -Holder continuous for any 7 < H — 1/4. 

Proof. Assume £1 < £2- We then have 

X(t 2 , x) - X(t u x) = A(t u t 2 , x) + B(t 1 ,t 2 , x), 



with 



and 



Hence 



A(ti,t 2 ,x)= f [ [G t2 s(x,y) - G tl - s (x,y)]B(ds,dy) 
Jo Jo 

B(t u t 2 ,x)=f [ G t2 - s (x,y)B(ds,dy). 
Jti Jo 



E [\X{t 2 ,x) -X(t u x)\ 2 ] < 2 (E [A 2 {t l ,t 2 ,x)]+ E[B 2 (t u t 2 ,x)]) . (15) 
We first note that (jOJ) and Lemma f2. II imply 
£[5 2 (£i,£ 2 , 



x 



duds\s - u\ 2H ~ 2 \ G t2 s(x,y)Gt 2 -u(x,y)dy 
-1 Jo 



<c I 2 ds{t 2 -s)- 1 / 2 I 2 \s-u\ 2H - 2 du 
Jt\ Jt\ 

< c(t 2 - t x ) 2H ^. (16) 
Now we will concentrate on the estimate on E[A 2 (t%, t 2 , x)]. By (j^J), we have 



E[A 2 {t u t 2 ,x)] =c H [ [ 

Jo Jo 



h rtl duds 



\s — u 



2~2H 



C x (s,u), (17) 



with C x (s, u) defined by 

C x (s, u)= [G t2 - S (x, y) - G tl _ s (x, y)] [G t2 _ u (x, y) - G tl ^ u (x, y)) dy. 



Thus, invoking Lemma \2.1\ we obtain that, for a given a < 1/2, 

(t — t ) 2a 

°^ U) ~ C ( tl -u)^(li-s)^ DxM ' 

where 

|1— a.\ r~< /„ „,\ s~~i „,\|1— o: 



D x (s,u)= \G ta - t (x,y)-G tl -s(x,y)\ a \G t2 - u (x, y) - G h - U (x, y)\ a dy. 
Jo 

It is then easily seen that D x (s, u) can be bounded by a sum of terms of the form 

F x (s,u)= f G 1 -_ a s (x,y)G 1 T Z a u (x,y)dy, 
Jo 

with cr, t G {ti, ^2}- This latter expression can be bounded in the following way: 

F x (s,u) < Qf 1 G 2 £; a) (x,y)dyY G 2 T { l: a) (x,y)dy 

c 



< 



(tl - s)V4-a/2( tl _ u y/4-a/2- 

We have thus obtained that 

Hl rtl duds 



E[A 2 (h,t 2 ,x)} < C (t 2 -t!) 2Q r [ 

Jo Jo 



\ S - u\ 2 - 2H {h - s) 1 /4+«(t 1 _ u )l/4+a • 

Now thanks the change of variable v = ^— -, the latter integral is finite whenever a < 

ti ' ° 

H — 1/4, which, together with ()15|) and (jl6|) . ends the proof. □ 

3 Ito's formula for the heat equation 

Let us turn to the main aim of this paper, namely the Ito-type formula for the process X 
introduced in (|T2"|) . The strategy of our computations can be briefly outlined as follows: 
first we will try to represent X as a convolution of a certain kernel M with respect to W, 
with reasonable bounds on M. Then we will be able to establish our Ito's formula for 
a smoothed version of X, involving a regularized kernel M £ for e > 0, by applying the 
usual Ito formula. Our main task will then be to study the limit of the quantities we will 
obtain as e — » 0. 

3.1 Differential of X 

Before getting a suitable expression for the differential of X, let us see how to represent 
this process as a convolution with respect to W. 
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3.1.1 Representation of X 

The expressions (jZJ) and (jHJ) lead to the following result (see ^\ 
Lemma 3.1. Let tp be a function in \Ht\- Then 



o jo 



<p(s,y)B(ds,dy) 



t r l 



JO 



[ K H,T 1 [o,t]<f)(u,y)W(du,dy), 



with 



[ K H? 1 [0A ( P\( u >y) = 1 ia,t]( u ) / (p(r,y)d r K H (r,u)dr. 



Remark 3.2. This result could also have been obtained by some heuristic arguments. 
Indeed, a formal way to write |3J) is to say that, for t > and y e [0, 1], the differential 
B(t,dy) is defined as 

B{t,dy) = [ K H {t,s)W{ds,dy). 
Jo 

Thus, if we differentiate formally this expression in time, since K H (t,t) = 0, we obtain 



d t B{t, dy) 



d t K H {t,s)W(ds,dy) 



uo 



dt. 



Since dtKn{t, s) is not a L 2 -function, the last equality has to be interpreted in the following 
way: if (p is a deterministic function, then 



^0 



<p(s,y)B(ds,dy) 



t r l 



<p(s,y) 



o </o 

t r l 



JO 



W(du, dy) 



d s K H (s, u)W(du, dy) 



(p(s,y)d s K(s,u)ds 



ds 



which recovers the result of Lemma \3.1\ 

We can now easily get the announced representation for X: 
Corollary 3.3. The solution X to Ml\) can be written as 



X(t,x) 



t r l 



^0 



M t;S (x,y)W(ds,dy), 



with 



M tiS (x,y)= / G t - u (x,y)d u K H (u,s)du. 



(19) 



Proof. The result is an immediate consequence of the proof of Proposition ^. 31 and Lemma 

ECU □ 
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3.1.2 Some bounds on M 

The kernel M will be algebraically useful in order to obtain our Ito's formula, and we will 
proceed to show now that it behaves similarly to the heat kernel G. To do so, let us first 
state the following technical lemma: 

Lemma 3.4. Let f be defined onO<r<t<T by 



f(r, t) = j\t - uy^iu - r)~ a exp **, 

for a constant k > 0, x G [0,2] and a G (0,1). Then, there exist some constants 
Ci, c 2 , c 3 , c 4 > such that 

f(r,t) < dCt-rr^-^exp (20) 

and 

d t f(r, t) < c 3 (t - r )-^W exp f-^) • (21) 

Proof. Recall that, in the remainder of the paper, k stands for a positive constant which 
can change from line to line. Notice also that ([201) is easy to see due to 

f(r, t) < exp (-jZ^j J\t ~ u)- 1/2 {u - r)' a dr. 

Now we will concentrate on ([21)1: let us perform the change of variable v = j^. This 
yields 

f(r, t) = (t- r)- {a - l,2) J\l ~ v)-^v- a exp (~ (1 _ ^ _ r) ) dv, 

and thus 

d t f(r,t) = giir,t) + g 2 (r,t), 

with 

9 

KX 



gi ( r ,t) = Kx 2 (t-rr^^ I (i- v )-3/V a exp(- ^ )dv 



and 



»(r,«> = g " <») («" 0^ /V - *r"V«exp (- (1 _g_ r) ) *■ 
Therefore, thanks to the fact that u t— > ue~ u is a bounded function on IR + , we have 
gi(r,t) < c(t - r)-( Q+1 / 2 ) J\l - v)-^v-«e W (- _^ —^ jdv 

KX X ' 



2{t — r) ) J 
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which is an estimate of the form (|2U|). Finally, it is easy to see that 



U-Ar.t) < r(f- 21 -XP ( "^TyJ jf (1 - vr^v^dv, 



which completes the proof. 



□ 



We are now ready to prove our bounds on M: 

Proposition 3.5. Let M be the kernel defined at Ufy . Then, for some strictly positive 
constants C5, Cq, C7, eg > ; we have 



M t>s (x,y) 
< c 5 (t-s) 



-(l-iO 



H ~ 1/2r ( <*{x-vY \ , /- c 6 (x + y-2) 2 
exp + exp 1 



t-s 



t-s 



and 



\d t M t)S {x,y)\ 

< C7 ( t _ s )-(2-H) ^ 



H-l/2 



exp 



c 8 (>-2/) 2 ^ , f c 8 (x + ?/-2) 2 
+ exp 



t — s J \ t — s 

Proof. First of all, we will use the decomposition (fTljl . which allows to write 

M M (x,y)= / G 1>t - u {x,y)d u K H (u, s)du+ j R t _ u (x,y)d u K H (u, s)du. 

J s J s 

Now the result is an immediate consequence of Lemma EOI applied to a < | — H, the only 
difference being the presence of the term (u/s) 11 ^ 1 ^ 2 , which can be bounded by (t/s) 11 ^ 1 ^ 2 
each time it appears. This yields the desired result. 

□ 



3.1.3 Differential of X 

With the representation (jl8j) in hand, we can now follow the heuristic steps in Remark 
13.21 in order to get a reasonable definition of the differential of X in time. That is, we can 
write formally that 



X(dt, x) 



t rl 



d t M t , s (x,y)W(ds,dy) 



dt, 



Uo Jo 

which means that if (p : [0, T] x [0, 1] — > R is a smooth enough function, we have 



(p(t, x)X(dt, x) 



(f{t,x) 



t r l 



JO 



d t M ti s(x,y)W(ds,dy) 



dt 



T pi 



JO 



W(ds, dy) 



<f(t,x)d t M ttS (x,y)dt 
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Note that this expression may not be convenient because it does not take advantage of 
the continuity of (p. But, by Proposition 13. 5[ we can write 

T pT 

ip(t, x)d t M tiS (x, y)dt = J (ip(t, x) - ip(s, x))d t M t>s (x, y)dt + ip(s, x)M T , s (x, y). 

Here again, we can formalize these heuristic considerations into the following: 

Definition 3.6. Let tp : Q x [0, T] x [0, 1] — > R be a measurable process. We say that p 
is integrable with respect to X if the mapping 

{s,y)^[M^ x ip]{s,y):= [ {<p{t,x) - <p(s,x)) d t M t;S {x,y)dt + ip{s,x)M T , s {x,y) (22) 

J s 

belongs to Dom{8), for almost all x G [0, 1]. In this case we set 

ip{t,x)X(dt,x)= [ [ [M*^}(s,y)W(ds,dy). 
Jo Jo 

Remark 3.7. Just like in the case of the fractional Brownian motion £7J/ or of the heat 
equation driven by the space-time white noise [6], one can show that J Q T <p(t,x)X(dt,x) 
can be interpreted as a divergence operator for the Wiener space defined by X. 

Remark 3.8. It is easy to see that Proposition ^. 51 implies that ip : [0, T] — > R is integrable 
with respect to X if it is (3 -Holder continuous in time with f3 > 1 — H . 

3.2 Regularized version of Ito's formula 

The representation (fTHj) of X also allows us to define a natural regularized version X e of 
X, depending on a parameter e > 0, such that t i— > X £ (t,x) will be a semi-martingale. 
Indeed, set, for e > 0, 

M^ a (x,y)= / G t _ u+£ (x,y)d u K H (u + e,s)du, 

J s 

and 



X £ (t,x)= / Ml s (x,y)W(ds,dy). (23) 
Jo Jo 

We will also need a regularized operator M^'* (see defined naturally by 

MJV] ( s > V) = J (f( r , x ) - x))d r M £ r s (x, y)dr + p{s, x)M^(x, y). 

Our strategy in order to get an ltd type formula for X will then be the following: 

1. Apply the usual ltd formula to the semi-martingale 1 1— > X £ (t,x). 

2. Rearrange terms in order to get an expression in terms of the operator M^'*. 
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3. Study the limit of the different terms obtained through Steps 1 and 2, as e — > 0. 

The current section will be devoted to the elaboration of Steps 1 and 2. 

Lemma 3.9. Let e > 0. Then, the process 1 1— > X £ (t,x) has bounded variations on [0,T], 
for all x G [0, 1]. 

Proof. The Fubini theorem for W and the semigroup property of G imply 

X £ (t,x) = / / G t -v+%(x,z) ( / / Gf E / 2 (-z, y)d u K H (u + e, s)W(ds, dy) ) ctedu, 
Jo Jo ' \Jo Jo / 

and notice that this integral is well-defined due to Kolmogorov's continuity theorem. 
Therefore, since t i— ► G t _ u+e /2(x, 2) is also a C 1 -function on [u, T], we obtain that X £ is 
differentiable with respect to t G [0,T], and 

d t X £ (t,x) 

= / / d t G t - u+ e(x, z) ( / / G £ / 2 (z,y)d u K H (u + e, s)W(ds,dy) ) cfedu 
Jo Jo \Jo Jo / 

+ J G £ / 2 (x,z) {^J J G £ / 2 (z,y)d t K H (t + e,s)W(ds,dy) S jdz, 

which is a continuous process on [0, T] x [0, 1], invoking Kolmogorov's continuity theorem 
again in a standard manner. 

□ 

An immediate consequence of the previous lemma is the following: 
Corollary 3.10. Let t G [0, T], x G [0, 1] and e > 0. Then, 

d t X £ (t,x) = QT d t Gt- u+£ (x,y)d u K H (u + e,s)dujW(ds,dy) 

+ 11 G £ (x,y)d t K H (t + e,s)W(ds,dy) 



J0 JO 
t rl 



JO 



d t M[ s (x,y)W(ds,dy). 



Proof. The result follows from Fubini's theorem for W and from the semigroup property 
oiG. 

□ 

Now we are ready to establish our regularized Ito's formula in order to carry out Steps 
1 and 2 of this section. 

Proposition 3.11. Let f be a regular function in C^(R) ; e > 0, and X s the process 
defined by \2ty . Then, for t G [0,T] and x G [0, 1], M^* f'(X e ) belongs to Dom(5) and 

f(X £ (t, x)) = f(0) + A li£ (t, x) + A a , e (t, x), 
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where 



Ai, e (t,x) 



t rl 



JO 



(M%f(X*)) (s,y)W(ds,dy) 



is defined as a Skorohod integral, and 



with 



A 2)£ (t,x)= f f"(X £ (s,x))K e , x (ds), 
Jo 

rs rv 2 , 

K £ ,x(s) = J dv 2 J dv 1 G2(s+e)-v 1 -v 2 (x,x)\H{2H - l)\Vx - V 2 \ 



2H-2 



-d, 



Vl+E 



K H {v\ + e, u)K H (v 2 + e, u)du 



"i 



-d V2 {K H {vi + £,Vi)K H (v 2 }■ 



(24) 



Proof. By Corollary we are able to apply the classical change of variable formula to 
obtain 



f(X%t,x)) = f(0) + / f\X%s,x)) 



s r l 



d s MUx,y)W(du,dy) 



ds. (25) 



Uo Jo 

Moreover, the derivative of f'(X £ (s,x)) in the Malliavin calculus sense is given by 
D v>z [f(X £ (s,x))\ = Ml v (x,z)f(X%s,x))l {v < s} . 

Since the last quantity is bounded by c £ v^~ H for e > 0, then invoking formula (|10j) for 
the Skorohod integral, we get 



's rl 



s rl 



JO 



d s M s (x,y)W(du,dy) 



o Jo 



f\X £ (s } x))d s MUx, y)W(du, dy) 



s rl 



JO 



+ f"(X £ (s,x)) / (d a Ml„{x,y))Ml„{x,y)dudy. (26) 



Denote for the moment the quantity J* (d s M e u (x,y)) M^ u (x,y)dudy by h x {s). Then, 
combining (|23J) and (J2SJ), proceeding as the beginning of Section and applying 

Fubini's theorem for the Skorohod integral, we have 

f(X%t,x)) = f(0) + A lte (t,x)+ [ f"(X e (s,x))h x (s)ds. (27) 

Jo 

We can find now a simpler expression for h x (s). Indeed, since M^ s (x,y) = 0, it is easily 
checked that 



hJs) = ~<9 



LJo Jo 



{Ml u (x,y)) 2 dudy 



(2* 
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Furthermore, the semigroup property for G yields 



i 

(Ml jU (x,y)) 2 dudy 
u J o 

s ps ps pi 

du dv l dv 2 dyG s+£ _ vl (x,y)G s+£ _ V2 (x,y)d Vl K H (vi + e,u) 

Ju Ju JO 

■ d V2 K H (v 2 + e,u) 



and this last expression is equal to 

ps ps ps 

2 du dvi dv 2 G 2 { s +e)-v 1 -v 2 {x,x) {d vl K H {v x + e,u)) d V2 K H (v 2 + e,u) 

J0 J u J v\ 



ps PV2 

2 dv 2 dviG 2 ( s+£) 
Jo Jo 



(x, x)[ I (d Vl K H (vi + e, u)) d V2 K H (v 2 + e, u)du 



(29) 



But 



I (d vi K H (vi + e,u)) d V2 K H (v 2 + £,u)du 
Jo 

I K H (v 1 + e,u)K H (v 2 + e,u)du 
Jo 



= d d 

-d V2 [K H (vi + e, vi)K H (v 2 + e, vx)] 

r pvi+e 

= H(2H -l)\ Vl -v 2 \ 2H - 2 -d V2 8 Vl / K H (v 1 + e,u)K H (v 2 + e,u)du 

.J V\ 

-d V2 [K H (vi + e, vi)K H (v 2 + e, Vi)) . 
By putting together (}2"9"j) and (J30J), we have thus obtained that 

1 f s f 1 

{M! Jx,y)) 2 dudy = K SiX {s), 



(30) 



o Jo 



where K EiX (s) is defined at (j2U). By plugging this equality into (|27jl and (|2"S)l. the proof 
is now complete. 

□ 



3.3 Ito's formula 

We are now ready to perform the limiting procedure which will allow to go from Proposi- 
tion E^2t the announced Ito formula. To this end we will need the following technical 
result, which states that the modulus of continuity of t i — > X £ (t,x) can be bounded from 
below by any v < H — 1/4, independently of e. 

Proposition 3.12. Let X e be given by (H3J). Then for t u t 2 e [0,T] and x e [0, 1], there 
is a positive constant c ( independent of e) such that 

E(\X s (t 2 ,x)-X e (t 1 ,x)\ 2 ) <c|t 2 -*i| 

for any v < H — -r. 



I 2;/ 
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Proof. Suppose that ti < t 2 . Then 



E(\X%t 2 ,x)-X%t l ,x)\ 2 ) < 2/ 1 / (MI 2iS {x,y)-Ml s {x,y)fdyd t 

Jo Jo 



'ti JO 

Now using the fact that d u K H (u, s) > 0, we have 
(M f e 2> ,(x,y)) dyds 

nl / Z^+e \ 2 

G t . 2 +2s-u(x,y)d u K H (u,s)duJ dyds 

-t2+e fl / ft2+£ 

< 



+2 / <2 f 1 (Ml s (x,y)) 2 dyds. (31) 




y fy l[t 1+£ ,t 2+e ](M)Gt 2+2e - u (a;,?/)9 u K H (M, s)c/mJ dyds 
= H(2H-1) / / |u - v\ 2H ~ 2 G t2+ 2e-u(x,y)G t2 +2e-v(x,y)dydudv 

Jti+e Jti+e Jo 

<c(t 2 -ti) 2ff -5, (32) 

where the last inequality follows as in (jlfij) . 

On the other hand, it is not difficult to see that 




./o 



< 



2/) - G tl y)} d u K H (u, s)du dyds 

\Js+e J 



I / G t2+ 2e-u(x,y)d u K H (u, s)du J dyds 

\Jh+s J 

B 1 + B 2 . (33) 



Observe now that we can proceed as in (j32|) to obtain 

B 2 <c(t 2 -t 1 ) 2H -*, (34) 

and it is also readily checked that 



y)|<9 n fT H (w,s)dw dyds 

= 2H(2H-1) [ 1+£ [ 1+£ [ \u-v\ 2H - 2 \G t2+2£ „ u (x,y)-G tl+2e „ u (x,y)\ 
Jo Jo Jo 



■\G t2+2E -. v (x,y) - G tl+2e -v(x,y)\dydudv. 
Finally, the proof follows combining (|17p. and ()31|)- 



□ 
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Let us state now the main result of this paper. 

Theorem 3.13. Let X be the process defined by and f e C 2 (R). Then, for t E [0, T] 
and x e [0, 1], the process M^ x f'(X) belongs to Dom(5) and 

f(X{t, x)) = /(0) + A ± {t, x) + A 2 {t, x), 

where ^ i 

A 1 (t,x)= [ [ (MU'(X))(s,y)W(ds,dy) 
Jo Jo 

and t 

A 2 (t,x) = ^J f"(X(s,x))K x (ds), 

with 

PS PS 

K x (s) = H(2H - 1) / / G 2s - Vl - V2 (x,x)\v 1 - v 2 \ 2H - 2 dv 1 dv 2 . 
Jo Jo 

As mentioned before, in order to prove this theorem, we use the regularized Ito formula 
of Proposition I3.11| and we only need to study the convergence of the terms Ai j£ and 
A 2j e appearing there. However, this analysis implies long and tedious calculations. This 
is why we have chosen to split the proof of our theorem into a series of lemmas which will 
be given in the next section. 



3.4 Proof of the main result 

The purpose of this section is to present some technical results whose combination provides 
us the proof of our Ito's formula given at Theorem 13.131 We begin with the convergence 
A2, e — > A 2 , for which we provide first a series of lemmas. 

Lemma 3.14. Let L\ be the function defined on [0, T] by 

rs pV2 

L\{s) = dv 2 dv 1 G2( s+ e)-vx-v 2 (x,x)K H (v 1 + e,vx)d V2 K H (v2 + £,v 1 ). 
Jo Jo 

Then s 1— > d a L\(s) converges to in -^([O, T]), as e j 0. 

Proof. Note that by (fPfj) we only need to study the convergence of d s L £ u (s), where 

L ii( s ) = [ dv 2 [ dvi + -_d V2 K H (v 2 + e,vi). (35) 

Jo Jo y/2(s + e) - v t - v 2 

Indeed, this term will show us the technique and the difficulties for the remaining terms. 

We will now proceed to a series of change of variables in order to get rid of the 
parameter s in the boundaries of the integrals defining L\ x : using first the change of 
variable z = V2 ~ Vl . and then 9 = Vi/s, we can write 



L £ n (s) = c H s 3 - 2H / (u-s9) H -lu H -tdu 



sO+e 



3 it 1 



'1 - 



Q2H-1 

f 1 (e + sO + zs(l - 6)) H -^ M m Nff3j f/1 
/ v - y >> ( Z s(l-9)+e) H -*dzd6. 

Jo y/2e + s(l-6)(2-z) 
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Hence, the change of variable v = u — s6 leads to 



02H-1 



f 1 (£ + s9 + zs(l -0))*"*, M n . . H _ 3j Jn 
/ , v " - (zs(l-6)+e) H idzdO 

Jo J2e + s(l-8)(2-z) 



/ 2e + s(l-9)(2- z) 

Therefore, by differentiating this expression in s, we end up with a sum of the type 

5 

where 



- 1 {e + s6 + zs{l -6)) H -^ 



1 'zs(l-6)+e) H ~ldzd6, 



y/2e + s(l - 9)(2 - z) 

and where the terms L\ l2 , . . . , L £ ll5 , whose exact calculation is left to the reader for sake 
of conciseness, are similar to L e m- 
Finally, we have 

Lln(s) < c H s~ H (jfV-**,) (T {l ~ 2 uT de ) f * H ~ ldz < c H e H - 1/2 s- H , 

and it is easily checked that this last term converges to in L 1 ([0,T]). Furthermore, it 
can also be proved that |-^fi,-(s)| < cL £ lu (s) for 2 < j < 5, which ends the proof. 

□ 

Lemma 3.15. Let L 2 be the junction defined on [0,T] by 

rs i-V2 / rvi+e \ 

L|(s) = / dv 2 / dviG 2 ( s+E )- Vl - V2 (x,x)d vl I / K H {v\ + e, u)d V2 K H (v 2 + e, u)du ] . 
Then s i— > 9 s L|(s) converges to in L 1 ([0,T]) ; as £ | 0. 

Proof. As in the proof of Lemma IH. 141 we only show the convergence of d s L 2l (s), where 

PS PV'2 J 

L 2i( s ) = / ^ 2 / ^i 7= . 9 w L(«i,t;2), 
Jo Jo v 2 (s + £J — v x — v 2 

with 

L(vi, v 2 ) = / K H (vi + e, u)d V2 K H (v 2 + e, u)du. 

J VI 
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Towards this end, we will proceed again to a series of changes of variables in order to 
eliminate the parameter s from the boundaries of the integrals: notice first that the 
definition of Kh, and the change of variables 9 = and z — r — v\ yield 



L(vi,v 2 ) =c H (v 2 + e) H -5 / (v! + z) H -^ / (V! + 6z) 



Thus 

d Vl L(v 1 ,v 2 ) = (v 2 + e) H ~^ 



A-2H 



z H -\ 



-H 



{y 2 + e-vi- 9z) H ~^d9dz. 



c H / Oi + z) H ~l / ( Vl + Bz) 



A-2H 



-H 



- c H I ( Vl + z) H -? I O + 9z) 



(«i + 9z) 



(l-9)i~ H 



• (v 2 + e - Vi - 9z) H '^d9dz 
■•- 2U Z " 2 (v 2 + s-vi -9z) H -'U9dz 



+ C H I (vi + z) H 3 



1-2H 



Z 2 



(1- 



-H 



(v 2 + e-vi- 9z) H ~^d9dz 



(36) 



Hence, it is easily seen that L\\ is a sum of terms of the form 



Qa,/3,v (s) 



V2 



dv 2 / dv\ 



(v 2 + e) H ^ 



a/2(s + e) - vi - v 2 Jo 



(v 1 + BzY- 



-H 



1 /■! 



s drj du 



o Jo 



77(577 + e) H 2 
a/2(s + e) - usrj - sr] Jo 



(usrj + 9zf- 



(vi + z) a 

(v 2 + e-vi- Qz) v dBdz 
(usrj + z) a 

(si] + e — sr/u + 6z) 1 'd9dz, 



-H 



by applying the changes of variable u = vi/v 2 , and rj = v 2 /s. Differentiating this last 
relation, we are now able to compute c^L^s), and see that this function goes to as 
e I in L 1 ([0,T]), similarly to what we did in the proof of Lemma f3. 141 

□ 

Lemma 3.16. Let L\ be the junction defined on [0, T] by 

PS PV2 

Ll(s)=H(2H-l) / dv 2 / dviG 2is+£ y Vl . V2 (x,x)(v 2 - Vl ) 2H - 2 . 
Jo Jo 

Then <9 s Lf(s) tends to \K x (ds) in L 1 ([0,T]) ; as e [ 0. 

Proof. As in the proofs of Lemmas 13.141 and 13.151 we only need to use the change of 
variables z = V\jv 2 and 9 = v 2 /s. 

□ 
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Lemma 3.17. Let X and X s be given in and \2'J\) . respectively. Then X £ (-,x) 
converges to X(-,x) in L 2 (Q x [0, T]) and, for t G [0,T], X e (t,x) goes to X(t,x) in 
L 2 (n), as e 10. 

Proof. The result is an immediate consequence of the definitions of the processes X e (-, x) 
and X(-, x), the fact that \Mf s (x, y)\ < c(t — s) H ~ 1 s2~ H and of the dominated convergence 
theorem. 

□ 

We are now ready to study the convergence of the term A 2e : 
Lemma 3.18. Let t G [0, T] and x G [0, 1]. Then the random variable 
B £ 2 (t,x) 

f"(X%s,x))d s U^ dv 2 J dv 1 G 2{s+E y Vl _ V2 (x,x)(v 2 -v 1 ) 2H - 2 )ds 

converges to A 2 (t,x) in L 2 (Q) as e j 0. 
Proof. Since f" is a bounded function, then 
E(\B e 2 (t,x)-A 2 (t,x)\ 2 ) 

< c f E {(f"(X(s,x)) - f"(X%s,x))) 2 ) \d s K x (s)\ds 

\d s K x (s) - H(H -^)d s dv 2 j dviG 2 ( s+£) - Vl - V2 (x, x)(v 2 - vi) 2H ~ 2 \d^j . 



+ c' 



Hence the result is a consequence of Lemmas 13.161 and 13.171 and the dominated conver- 
gence theorem. □ 

Now we study the convergence of A l e to A x in L 2 (Q). 

Lemma 3.19. Let X and X s be given in HT^) and p?3J) . respectively. Then, for t G [0, T] 
and x G [0, 1], 

E Qf jf [(MU'(X))(s,y) - (M E t ;:nX £ ))(s,y)] 2 dyds^ -> 

as e [0. 

Proof. We first note that 

E Qf J [(MU'(X))(s,y) - (M° : :f(Xe))(s,y)] 2 dyds 
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can be bounded from above by: 
ft pi r pt 



cE 



jo 



(f(X(r,x))-f(X(s,x)) 



f (X £ (r, x)) + f'(X%s, x)))d r M r>s {x, y)dr 



dyds 



+ cE 



+ cE 



+ cE 



t pi 



Jo 
t pi 



(f\X%r,x))-f\X%s,x)))(d r M r , s (x,y)-d r MUx,y))dr 



dyds 



[(f(X(s, x)) - f'(X £ (s, x)))M tiS (x, y)f dyds 



JO 

t pi 



JO 



[f(X B (s, x))(M M (x, y) - Ml s (x, y))} 2 dydt 



c{B x +... + B 4 



Next observe that 



(37) 



Bo < 



t pi 



o Jo 



E ((f(X%r,x)) - f'(X%s,x))) 2 ) \d r M r , s (x,y) - d r M^ s (x,y)\dr 



dyds. 



\d e M 0!S (x,y)-d Ml s (x,y)\de 

Now notice that Proposition 13.121 and the inequality 

\d r M*{x,y)\ 



< c 



r + e 



(r — s + e 



<H-2 



exp —C\ 



(g - y) 2 

s + (r — s) 



+ exp -ci 



{x + y-2f 



e + (r — s) 



imply, for (3 small enough, that 

E ((f'(X%r,x)) - f'(X%s,x))) 2 ) \d r M r , s (x,y) - d r M^ s (x,y)\ 

goes to as e | and that it is bounded by cs^~ H {r — s) 311 ^^^ 13 . Thus 

B 2 ^0 

because of the dominated convergence theorem. 
Since /' is a bounded function, then 



(38) 



Ba < c 



t pi 



JO 



(M t>s {x,y)-Ml s {x,y)) 2 dyds, 



which goes to due to the definition of M £ and the dominated convergence theorem. 
Hence, by (J37J), and (|3*H|). we only need to show that B\ + B 3 — > as £ I to finish the 
proof. This can been seen using Lemma 13.171 and proceeding as the beginning of this 
proof. □ 
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Lemma 3.20. Let X and X s be given by KTty) and W^l . respectively. Then, fort E [0, T] 
and x G [0, 1], M^ x f'(X) belongs to Dom (5). Moreover 

5 {M%f{X*)) - 5 (MtJ'(X)) 

as e 10 in L 2 (tt). 

Proof. The result follows from Lemmas l3.14H3.lSl and from the fact that 5 is a closed 
operator. □ 
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